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A new super symmetric approach to the analysis of dynamical symmetries 
for matrix quantum systems is presented. Contrary to standard one dimen- 
sional quantum mechanics where there is no role for an additional symmetry 
due to nondegeneracy, matrix hamiltonians allow for non-trivial residual 
symmetries. This approach is based on a generalization of the intertwining 
relations familiar in SUSY Quantum Mechanics. The corresponding matrix 
supercharges, of first or of second order in derivatives, lead to an algebra 
which incorporates an additional block diagonal differential matrix operator 
(referred to as a "hidden" symmetry operator) found to commute with the 
superhamiltonian. We discuss some physical interpretations of such dynam- 
ical systems in terms of spin 1/2 particle in a magnetic field or in terms of 
coupled channel problem. Particular attention is paid to the case of trans- 
parent matrix potentials. 
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1. Introduction 

Supersymmetric Quantum Mechanics (SUSY QM) is an interesting framework 
to analyze non-relativistic quantum problems. In particular it allows to investigate 
the spectral properties of certain quantum models as well as to generate new systems 
with given spectral characteristics. 

In general it is well known that SUSY algebra provides the relation between (su- 
per)partner hamiltonians (which are often referred as "bosonic" and "fermionic") 
associated to a two-fold degeneracy of levels. 

Much less attention has been paid in the literature to the possibility of using 
SUSY as a tool to study individual ("internal") symmetries of each superpartner 
hamiltonian. In principle this can be useful if it allows to integrate partially dynamical 
systems by discovering an additional dynamical symmetry which we shall refer to as 
a "hidden" symmetry. We will show that indeed hidden symmetries may be revealed 
by SUSY inspired approaches 0. 

The standard SUSY QM relations read 

{Q\Q-} = H=( h ^ h % ] = ( *r? ); (i) 

h® = -d 2 + V [i \x)- q ± ~- 

Q-=(°-l) Q + =ln\)> (2) 

{Q ± ,Q ± } = o. (3) 

[£±#1 = 0; h W q + = q + h {2) (4) 

In general there are different realizations of this algebra, e.g. multidimensional 
and matrix ones || . It is also possible to generalize the algebra itself by preserving 
Eqs. (0), (f|) and allowing, by a non standard form of the intertwining operators q ± , 
to modify (|TJ) to become 

where the diagonal operator K is no more the superhamiltonian but has in general 
the nature of a symmetry operator 

[k {i \h {i) ] =0. (6) 

^Notice that the existence of a symmetry operator is not necessarily associated with degeneracy of 
eigenstates of the Hamiltonian. One can convince oneself of this fact for matrix Quantum Mechanics 
taking 2x2 diagonal hamiltonian with components having different spectra and no degeneracy. In 
this case a symmetry operator is obviously 03 playing the role of a grading operator. 
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This generalization was already discussed for 1-dim ||] and 2-dim H Quantum 
Mechanics. For 1-dim systems and intertwining operators of second order in deriva- 
tives the only relevant case was the one for which K is a function of H, 

K = H 2 — 2aH + (3 

where a and (3 are constants. For the 2-dim case there exists the possibility of having 
a central charge R, which commutes with all elements of algebra, such that 

K = f(H) + R. (7) 

A new, supersymmetrical, method was elaborated ||], || to investigate hidden dy- 
namical symmetries of quantum systems. The existence of such a differential operator 
R implies a dynamical symmetry (unknown a priori) made apparent by the general- 
ized SUSY algebra. 

Let us remind the physical impact of the supersymmetric approach to the Dirac 
equation || with some applications to superconductivity Jj]], to pseudorelativistic be- 
haviour of electrons in two-band systems || and to attempts to a diagonalization 
procedure 0) in nuclear and atomic physics coupled channel problems and finally in 
the treatment of particles with spin in an external magnetic fields |§. Therefore it is 
important to discuss the role of symmetry operators for 1-dim dynamical systems of 
Schrodinger type with matrix potentials. Formally the most straightforward method 
is to study the commutator of the Hamiltonian matrix with a generic differential op- 
erator matrix R and solve the corresponding system of differential equations obtained 
by imposing the commutator to vanish [h, R] = 0. Incidentally one can easily check 
that for the scalar case this system of equations (after possible subtraction of the 
Hamiltonian) does not allow for non-trivial solution. On the contrary in the matrix 
case non-trivial solutions exist even for the case of symmetry operators of first order 
in derivatives. For higher order derivatives the equations become rather cumbersome 
and it is too difficult to provide the general discussion of the solutions. 

For this higher order case a method of solution suggested by supersymmetry seems 
to be useful. It starts from the same idea of factorization originally proposed by 
Schrodinger [|llj with the related intertwining operators (see ((D) but now it is not 
applied to the Hamiltonian but rather to the symmetry operator: this method reduces 
the order of differential equations one must solve without however increasing their 
number. In this paper we shall find non-trivial genuine second order operators R even 
for 1-dim first order intertwining matrix operators. 

The paper is organized as follows. In Section 2 we shall investigate the SUSY 
approach with intertwining operators of first order in derivatives. A variety of matrix 
systems allowing for the genuine second order symmetry operators will be obtained. 
In the Section 3 this method will be applied to a more complex case of higher order 
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intertwining operators, examining in particular the conditions for the factorizability of 
these intertwining operators (reducibility) . The type of physical systems which we can 
describe in our formalism include a particle with magnetic moment in a magnetic field 
and more generally a class of systems with one continuous and one discrete degrees 
of freedom. A novel construction is given for transparent matrix potentials which are 
not duplications of standard scalar transparent potentials and are not generated by 
iterations of first order Darboux transformations. 

2. First order Matrix SUSY Quantum Mechanics 

We start from the general first order (in derivatives) representation of the com- 
ponents of supercharges in the case of 1-dim Quantum Mechanics 

q + = Ad+B; q- = -A ] d + B\ (8) 

where A and B are matrices. Imposing that these operators intertwine the Hamilto- 
nians h^ 1 ' and (see Eq.(fD) reads 

(-d 2 + V w (x))(Ad + B) = (AO + B)(-d 2 + V {2) (x)), (9) 

where V {l \V {2) are hermitean potential matrices. Eq.@ amounts to solve the fol- 
lowing three equations: 

A' = 0; (10) 

V {1) A - AV {2) = 2B'; (11) 

- B" + V {l) B - BV {2) - AV (2) ' = 0. (12) 

The first equation implies that from now on we will assume A to be a constant matrix. 
As usual in the frameworks of SUSY QM the intertwining relations Eq. ([J) lead to the 
connection between the column eigenfunctions of Hamiltonians and : 

^ (1) (x) = (A3 + B)¥ 2 \x); ¥ 2 \x) = (Ad + B)^ {1 \x) (13) 

The Eq.([llJ) allows sometimes for zero-modes, and the spectra of the partner hamil- 
tonian coincide up to these zero-modes. 

What is not standard is the fact that Eq.(^) does not hold because the products 
of the supercharge components are no more equal to the Hamiltonians: 

q + q - = k w = -AA^d 2 + (AB ] - BA ] )d + BB ] + AB V ; (14) 



q-q+ = k {2) = -A ] Ad 2 + (B ] A - A ] B)d + B ] B - A ] B' '. 



(15) 
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It is natural to consider separately det A ^ and det A = because the way to 
solve the system of equations (|T0|) , ([TlD , ([12]) differs for the two cases. In the case 
det A ^ it is possible to subtract the Hamiltonians and from k^ 1 ' and k^ 2 \ 
being left with a symmetry operators of first order in derivatives; otherwise we deal 
with second order operators. Furthermore by a suitable similarity transformation 
induced by the matrix A itself it is possible to consider A = — 1 as a representative of 

the case det A ^ whereas we can take the matrix A — ( ^ g ) as representative 



of det A = 0. 

2.1. The case A = — 1, symmetry operators of first order 
In this case Eqs.(|TT|) and ( |1"2"D become 



y(2) _ y(l) = 2B >. ( 16 ) 

- B" + V {1) B - BV {2) + 1/ {2) ' = 0. (17) 
One can parametrize the potential matrices taking into account flTGp: 

y(l) = ^2 _ + ^y. y (2) = ^2 + ^/ + ( lg ) 

Then Eq . C|X7|) becomes simply 

W'(x) = [B(x),W(x)}. (19) 

Whoever is familiar with ordered products, e.g. in Gauge Field Theory, will not have 
difficulty to recognize that Eq.(|I§|) has not a simple local solution unless 

[B(x),B(y)]=0. 

Since effectively this last condition reduces the matrix problem to a scalar problem 
we do not find this case of any interest for us. 

^From (|16D one deduces that the non-hermitean part of B does not depend on x 
because of the hermiticity of the potentials. Therefore we can parametrize B(x) = 
B(x) + iC, where B and C are both hermitean matrices and C is a constant matrix. 
Correspondingly one has to solve the system of matrix equations 

W = [B, W] + i[C, W}; (20) 

W -W ] = -2i{B,C} (21) 
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which are the consequence of ([D]) and of the hermiticity of both potentials. Restrict- 
ing to 2 x 2 matrix potentials one can study the system of equations ( p0|) and (|2l|) 
by expanding all matrices in Pauli matrices and unity: 

B = b + bi<Ji; C = c + Ci<Ji; W = w + Wi(Ji, (22) 

with components bo, bi and cq, q real. The related symmetry operators read: 

RW = q + q - - h {1) = liCd + 2C 2 - 2iBC - W; (23) 

#( 2 ) = q -q+ - h {2) = 2iCd + 2C 2 - 2iCB + W. (24) 

While it is not restrictive to set c\ = c 2 = 0, it is not interesting to choose also 
C3 and Co both to vanish. Indeed in this case the symmetry operators above are no 
more differential operators and become proportional to a constant matrix W. 

A general solution of the nonlinear system of matrix equations (^0|), ([H]) amounts 
first to find a solution of the subsystem: 

Imw = -2(b c + b 3 c 3 ); 
Imw\ = —2bic ; 
Imw 2 = -2b 2 c ; 
Imw 3 = -2(6 3 c + 60C3); 

w = const; 
Rew 3 = const, 

however a complete solution of (|20|) , (f2l|) cannot be written. Particular solutions of 
( |2"DD and ( |2"1"D in terms of the components ( P2"D can be obtained by making specific 
ansatze: 

1) Let us assume 

c o — ^3 — Rew 3 = 0; Imw 3 (x) = 0; c 3 ^ 0, 

b 3 7^ would correspond to a trivial solution W(x) = Const. Then one obtains 

6 = 0, Rewi = /?cos(2c 3 x), Rcw 2 = — /9sin(2c 3 x), b 2 (x) = — bi(x) tan(2c 3 x), 

where (3 is a constant parameter and b\{x) = b\(x) cos(2c 3 a;) with b\{x) being an 
arbitrary nonsingular function. Correspondingly the potentials read 

V m2) (x) = bl(x)-cl + Rew 

+ [p cos(2c 3 x) =F b[(x) cos(2c 3 x) ± 2c 3 6i sin(2c 3 x)] ■ <j\ 

— [(3 sin(2c 3 x) =F b[( x ) sin(2c 3 x) =F 2c 3 b\{x) cos(2c 3 x)] • a 2 (25) 
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and the symmetry operators (|23|), ( |24| ) after subtraction of constants become: 

^(i),(2) _ iu^Qzs^^h^x) sin(2c 3 x)H — — cos(2c 3 x)]-ax^\bi{x) cos(2c 3 x) — — sin(2c 3 :r)]-a"2- 

2c 3 2c 3 

(26) 

Among the different interpretations concerning the physics of the matrix potentials 
one consists of a spin 1/2 neutral particle in a (inhomogeneous) magnetic field. 
It is necessary to assume that the magnetic field depends only on the coordinate 
x = x 3 and lies in the (x\, x 2 ) plane in order to ensure the vanishing of its divergence 
diBi = 0. While the motion in the (x±, x 2 ) plane is trivially free the dynamics is still 
rather interesting because of the x 3 motion |{J. Physicswise the inhomogeneity of the 
magnetic field is determined in this case by the requirement of vanishing of the scalar 
potential in (|25f) (neutrality of the particle). 

2) Let us assume 

c o — ^3 — Rew 3 = 0; Imw 3 (x) 7^ c 3 7^ 0, 

as before 63 7^ would correspond to a trivial solution W(x) = Const. Then one 
obtains 

1 (Rew 2 )' + 2c 3 Rew 1 -{Rew x )' + 2c 3 Rew 2 

b (x) = -—Imw 3 (x)] bi(x) = — ; b 2 (x) 



where 



2c 3 ! 2Imw 3 ' 2Imw 3 

Rewi(x) = J (Imw 3 ) 2 + (3 cos f(x) = \J (26 c 3 ) 2 + (3 cos f(x), 



Rew 2 (x) = J (Imw 3 ) 2 + (3 sin f(x) = \J (26 c 3 ) 2 + (3 sin f(x). 

(3 is a constant parameter as well as c 3 and f(x) and Imw 3 (x) are arbitrary functions. 
The potentials and the symmetry operators read: 

V^ 2 \x) = b 2 Tb' + bl + bj-cl + Rew 

+ [2Mi + Rewi =F b[] ■ (Ji + [2b b 2 + Rew 2 T b 2 ) ■ a 2 (27) 

R^ = ia 3 d ± \b 2 {x) - ■ a, ± [-b x {x) - ^] • a 2 . (28) 

zc 3 zc 3 

In terms of the "magnetic" interpretation given above one can now notice that the 
absence of the scalar potential in ( |2"TD is less restrictive because the magnetic field still 
depends on one arbitrary function. The intrinsic "periodicity" of the magnetic field 
forces a similar periodicity of the wave function. We warn however not to interpret the 
periodicity too naively since it depends in general on the properties of the arbitrary 
function f(x), e.g. asymptotically constant magnetic field can be incorporated in this 
scheme. 
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One can also find solutions for other ansatze like e.g. 

c = 0, c 3 7^ 0, Rew 3 = const ^ 0, 

similarly to 1), 2). 

2.2. The case detA = 0, the symmetry operators of second order 
The constant (see (|10D) matrix A has now the form 

Hoo) (29) 

We write the matrices V^'^ and B explicitly and for simplicity we assume them to 
be real: 

V- = ( S 11 ) (30, 



1 




(0 


(* 


bi 


b 2 


h 


by 



Equations (|TT| ) and ([12]) can be solved and the potentials can be written in a 
similar way as before in terms of Pauli matrices, however the reality condition forces 
the absence of a 02 term R. The hidden symmetry operators are now of second order 
in derivatives. 

In the case b 7^ it is possible to find several solutions dependent on arbitrary 
functions; in the case b = 0, a = 1 the form of the solutions simplifies. One solution 
is given by: 

v (1 \x) = -2cb' 2 b 2 2 ; 

bi ,2 ,2 c 2 ,/ 
— + &I + &2 + 72 -b[ + cx; 
1 b 2 

v®{x) = -2b' 2 ; (32) 

c 2 

where &i(x), bi{x) are arbitrary functions, 63 (x) = c-b 2 1 (x) and 64 = 0. The symmetry 
operators iJW'v 2 ) can be straightforwardly derived according to (p!4[), (p~5| ) and are 
second order differential operators. 

e The interpretation of the dynamics as a magnetic interaction of s = 1/2 particle is now still 
possible but we can also consider it as a coupled channel problem as discussed in S . 



v[ 1] 




= K + bl + c 7 


6 2 2 + ci; 






- -0+ 


/&'\ 2 5 

J + " 




(x) 


= -&; + ^ + 




4 2) 


(x) 


62 &2 
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For the second solution Vi'^ 2 \x) and v^^ 2 \x) are the same as in fl32|) and 



" — r °i + °fc + c i> 



" 2 W = 6s 



fe=l 
4 



t ,f W = M + ^l + ^ + 6 ; + I: ^ + c - 1 , 

where 62(2^), &3(^) are arbitrary functions such that 62 • &3 7^ 0, 64 = const 7^ and 



&i(*) = (Ms)" 1 
with a a constant parameter. 



26 4 



-(26 4 )- i (^ + 6^) + a 



3. Second order Matrix SUSY Quantum Mechanics 

Let us define the second order differential operators 

q + = ( g -)t = d 2 -2F(x)d + B(x). (33) 
q- = ( g +)t = 52 + 2^(3)0 + S t( x ) + 2F' t (a;). (34) 

where F(x) and -B(x) are 2x2 matrices. This representation for q can be inserted 
into (D-d). 

The intertwining relations are equivalent to a system of three non-linear matrix 



differential equations: 

y{l) _ y(2) + 4F / = (35) 

F" - + FV (2) - 5' - V (2) ' = (36) 

B" + V (2) " - V W B + - 2Fy( 2 )' = (37) 



Our attitude towards the solution of this system of equations is that we consider 
q ± , hS 2 \ h^ 1 ' to be essentially unknown except for Schrddinger form of hamiltonians 
and assumption of structure (|33|) of the supercharges, so the problem is to find the 
solution in terms of the matrices F(x), B(x), V^(x). 

Due to the complexity of the problem ( matrix, 2nd derivatives, non-linearity ) it 
does not seem realistic to search for a general solution in analytic form, instead we 
believe that techniques of Higher Order SUSY QM as developed in [§J can provide a 
useful tool for solving in an "indirect" way by the ansatz of factorizability of q^, i.e. 
restricting to the reducible matrix Higher Order SUSY QM. 
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Another possibility which we mention is a particular solution for which the terms 
appearing in Eq . (|36|) reduce to {F, F'} with the aim of a "direct" 
integration of this equation. A sufficient condition which allows this integration is 

y(l) + y(2) = 2 p( F ) 

where P is an arbitrary "scalar" function of the matrix F(x) like for example P = 
T>c n (x)F n (x) where scalar functions. 



3.1. Reducible Higher Order Matrix SUSY QM 

A specific ansatz consists in the factorizability of the operators of (|33"D, (53) in 
terms of ordinary superpotentials W(x) and W(x) : 

q + = q + q + = (-0 + W(x)) (-0 + W(x)), (38) 

connected by the ladder equation 

q-q+ = q +q 2 + A or W + W 2 = -W' + W 2 + A, (39) 

with A a constant hermitean matrix as will be clear later on. Let us assume fur- 
thermore A to be diagonal^. Then F and B of Eqs. (|33D , (|34[) are determined by the 
superpotentials W(x),W(x) : 

2F = W + W, B = WW - W' (40) 



The factorization Eq. fl3"8]) arises from two successive standard SUSY QM trans 
formations 

/ \ / q+qT \ 

h { q^qf I 



and 

( h \ = ( q+q^ + A 
^ h® + A J { g 2 -g 2 + + A 

by deleting the "intermediate" hamiltonian : 



(41) 
(42) 



H = f + A ) • (43) 

V q 2 q+ + A J 1 ; 

The matrix A has to be such that [H, = and consequently [qf, A] = 
which makes clear the reason why A has to be constant and that [W, A] = 0. We 

' If it is not it can be diagonalized by a constant unitary transformation which also affects other 
operators. 



Matrix Hamiltonians. 



12 



therefore can conclude that [h^ + A, A] =0 and this allows us to identify = A 
as a symmetry operator for h^ 2 \ If it further commutes with hP~> this operator is 
such that 

{Q + ,Q-} = (H) 2 -T-H (44) 
where T is a block diagonal matrix 

r -(ti)- w 

and thus it corresponds to a rather "trivial" operator. We shall from now on 
exclude such a case imposing that the operator 

R {1) = -qfAq^ (46) 

be non-trivial f\ The last case is rather interesting because it incorporates the possibil- 
ity of having genuine partner symmetry operators of different orders in the derivatives. 

In order to derive a solution we expand the previous operators in terms of Pauli 
matrices: 

W(x) =w + ww, W(x) =w + w 3 a 3 ; A = 5 + 5 3 a 3 ; 5 3 ^ 0. (47) 
To illustrate the techniques involved we give now an example: 

Example 1. 5 = 0; w (x) = Wo(x); w 2 (x) = 
Eq.(|39|) takes then the form 

2wqw 3 + w 3 — 2wqw 3 + w 3 — 5 3 = 

2w w 1 + w[ = (48) 

W 3 + w\ + — ^3 = 

with the following solution 

w = w = l/{2x), w 1 = l/x, w 3 = w 3 = 5 3 x/2. (49) 
In order to arrive to an interpretation of these results we have to write the potentials: 
V m (x) = W 2 -W = 7 /Ax 2 + 5 2 3 x 2 /A + (2/x 2 ) ■ a x ; (50) 
V {2 \x) = W 2 + W + A = -l/(4x 2 ) + 5 2 3 x 2 /A + 25 3 ■ a 3 . (51) 



9 In the concluding Section of H Eq. (mf) first appeared 
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e.g. eigensolutions These potentials contain centrifugal singularities and therefore 
we pose the eigenvalue problem on the semiaxis x > 0. The physical solutions 
L2— normalizable at the origin have the following behaviour: 

*<"(*) ~. + *«(*)- fi'fSV (52) 



*o y aix 5 / 2 — a 2 x 1 / 2 

where we display explicitly the diagonalization of matrix potential V^(x). Both 
potentials for £3 7^ lead to a discrete spectrum. 

The symmetry operator can be calculated from Eq. (fl6| ) and found to contain 
also o~ 2 type of terms: it can be written as 

R0) = a3 d 2 + (2i/x)a 2 d + {-8 3 ■ a x - {1/x 2 ) ■ a 2 + [l/(4x 2 ) - 5 2 x 2 /4] ■ a 3 }. (53) 

This operator also contains centrifugal singularities but one can easily prove that it 
maps L 2 — normalizable solutions ( |52|) into themselves. Therefore it is a true symmetry 
operator. The partner symmetry operator R^ = 63(73 is regular. 

As in Subsection 2.1. also this example allows an interpretation in terms of the 
external field as a magnetic field. It is possible to consider a magnetic field as a 
(pseudo)vector in a plane orthogonal to the one dimensional axis in which the particle 
is allowed to move x = x 2 []. Contrary to Subsection 2.1. the magnetic field in 
is homogeneous along the X3 axis, while in V^ 1 ) it is not homogeneous (depends on 
x 2 ) and has non zero components in the (2^, £3) plane. 

To describe other examples it is useful to introduce 

W(x) + W(x) = 2F(x) = 2/o + 2f,a t ; Ifx = w u 2/ 2 = w 2 
^9f) becomes: 

IF' - AF 2 + 2{F, W} = A. (54) 

The problem is now expressed in terms of a matrix non-linear differential equation 
( |54"D for F, W. Since we are unable to present a general (analytic) discussion we 
provide particular solutions of the system ([54]) which in terms of components can be 
rewritten: 

2/o - 4(/ ) 2 + Khf + 4(/ 2 ) 2 - 4(/ 3 ) 2 + 4/o^o + 4/3^3 = 5o; (55) 

f[ + 2w /i = 0; (56) 
f 2 + 2w f 2 = 0; (57) 



h This choice will allow to implement the condition of absence of sources of magnetic field auto- 
matically, diBi = 0. 
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2/s - 8/0/3 + 4/o^s + 4^0/3 = S 3 . (58) 
The solutions presented in the following will allow to construct explicitly the partner 
potentials and the symmetry operators by making use of the general expressions in 
terms only of the fa, (k = 0, 1, 2, 3) and Wq, w 3 : 

+ (2w w 3 - w 3 ) ■ a 3 ; 
V^(x) = w l + w l-w' + A{fl + fl + f Q )-5 

+ [2w w 3 -w' 3 + 8(2/0/3 - f w 3 - f 3 w ) + S 3 ] ■ 03. 

It is advantageous to use not the expression ( J4"6| ) directly but to define a symmetry 
operator by suitable subtraction and rescaling (see Section 1): 

~ (1) + 5 hW 0,,.,. r ^ 

R { ' = 7 = -(T3O + Al(f 2 (Tl - f 102)0 

^3 

+ (2w w 3 - w' 3 ) + 4(w 3< /i - iwofz) ■ <J\ + A(w 3 f 2 + iwofi) ■ o 2 
+ (w% + wl - Afl - 4/| - w' ) ■ 0-3, 

and we remind that = A. 

We now list two particular cases: 

Example 2. /„ = 0; f 3 £ 
The solution is of the type 

2/i (a) = 7iexp(-2 / w dx); 



2/2^) = 72exp(-2 / w dx); 

2/3(0:) = exp(-2 J w dx)[y 3 + <5 3 exp(+2 J w dx)}; 
w 3 {x) = (4/ 3 )-%-(7i 2 + 72)exp(+4 [ w dx) + Aft] 



Example 3. / 3 = 0; / ^ S = 
The solution of the Eqs.(|55"l) — (|58|) can be written as 
2/i (x) = 7iexp(-2 J w dx); 

2/2(2;) = 7 2 exp(-2 J w dx); 
w 3 (x) = $3/(4/o); 

fo(x) = l/2 v / 7l 2 + 7 2 2 exp(-2 J w dx) tanh(-y / 7l 2 + 7 | / exp(-2 / w rfx) +C). 
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3.2. Transparent matrix potentials 

We now explore a physical case for which the Hamiltonian describes free 
motion; in Eqs. (|35[) , ([36]), (|37D we assume = and, as a consequence of SUSY, 
becomes a so called transparent matrix potential [|j], 0: 

y(l) + AF > = o (59) 

F" - yWf - B' = (60) 

5" - = (61) 

These relations allow for further simplifications by eliminating and B in terms 
of F using Eq. (|59|) and differentiating Eq.(^) to get B" which is then introduced in 
Eq . (|6l|) leading to a linear algebraic equation for B |. Inserting back the expression 
for the derivative of B into Eq.(^) we obtain 

F" + AF'F + [(AF')-\F'" + AF' 2 + AF"F)\ = (62) 

or equivalently 

(F" + AF'F)' + 4F' /" (F" + AF'F)dx = (63) 

a non-linear 4-th order equation for F whose solution may allow to identify a class of 
transparent potentials in a 2 x 2 coupled channel problem. 

A sufficient condition for F to satisfy Eq. ([62]) is given by the simpler equation of 
second order 

(F" + AF'F) = jF' (64) 

with 7 an arbitrary constant number. It easy to verify the property that if F(x) is a 
solution of the nonlinear matrix equation (|64"D than F(x) = F(x)+a-I with a = const 
is again a solution of the same equation but for the shifted value of 7 = 7 + Aa. This 
peculiar property allows therefore to restrict to the case 7 = in (|64[) . The solution 
of this equation because of (|59D can be searched for, parametrizing F(x) = G(x) +iC 
with G(x) and C hermitean and C a constant matrix. We have thus to solve 

G"(x) + AG\x)G{x) + AiG\x)C = 0. (65) 

We can expand G and C in Pauli matrices, by a suitable rotation we can choose 
c i — c 2 — and furthermore we assume Cq = 0. Then we arrive to a system of 

1 Incidentally we note that Eq.([n]) is a Schrodinger equation for B with eigenvalue zero. This is 
related to the fact that — and therefore there exists a trivial zero-energy solution, namely the 
constant wave function. Letting A + act on this wave function one obtains the corresponding one for 
with the same zero energy as proportional to the matrix B(x) acting on a constant vector. 
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equations: 



g' + 2g 2 + 2g 2 = 2 7o ; 

91 + ^9o9i - 4£ij3C39j = 2 7i; i, j ■ 

eijk9j9k + 9oC 3 5 3i = 0; i,j,k 

93 = 273, 

with constant 7's and obvious meaning of (?'s and c's. 
A solution can be found for all 7^ = : 

1 



1,2; 
1,2,3; 



9o(x) = 

01 0*0 = 



4x + 

5-0 (x) cos0(:r); 
^(a;) = #o(a0 sin0(x); 
4>(x) = -4c 3 x + C, 



with (5 and £ arbitrary real constants. 

By shifting x — > rr — /3/4 and changing ( — > £ + 03/? we obtain the transparent 
hermitean matrix potential V^(x) with a singularity at the origin and asymptotic 
long range behaviour: 



V {1 \x) 



where 



x/ 



+ 



coscf)(x) — sin0(:r) 



r 



x 



•<7i + 



— - sin Mx) H — — cos <j>(x) 
x 2 x 



1 ~ A-ir-i 

; (l + N(x))-^a 3 N(x), 



ar 



X 



■ (7 2 

(66) 



JV(a;) = <7„ cos(4c 3 a;) + io 3 a n sin(4c 3 x) 
a n = U\ cos ( + <7 2 sin £. 
Due to the centrifugal singularity in this potential 

~ o ^(1 + + 0(1) 

one should consider the scattering problem on the semiaxis x > 0. 

To analyze the behaviour of wave functions at the origin it is useful to choose 
as a basis the following matrices: <J n ,a n = (02 cos £ — o"isin£ ) and a 3 with non- 

/ 1 \ _ / 1 \ / -i' 



standard realization: er„ 



-1 



1 



,<*3 







Then po- 



tential V^\x) becomes diagonal at the origin and its wave functions for x — > are: 



a • x 
b ■ x 
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The dynamical symmetry operator for hP> can be derived from 

R 1] = q + q- - {h w f = -AiCd 3 + {-AG 2 - AC 2 - Ai[C, G] - 2G') ■ d 2 - 
(8GG' + 8tCG' + 2G")d - (2G" { + AGG" + AiCG" + 16G' 2 ). 

This operator is also singular at the origin: 

R l) Aic 3 a 3 d 3 + (— a n - Ac\ + ldc 2 3 a n - 32eija;& n ) • d 2 + (^r<7 3 - 

x-*Q x X z 

Ac 3 , . n , 4ic 3 4c 3 12c 2 12c 2 32dj , . 

Nevertheless it maps wave functions regular at the origin into regular ones similarly 
to the Example 1 (Section 3.1.). 

Proceeding in the same way for 70 > 0,7^ = (k = 1,2,3) we find another 
solution: 

uj coth(4cjx + 0)\ 

iO COS0 

sinh(4u;x + 0) ' 

uj sin <fi 
sinh(4u;x + 0) ' 

with (f)(x) as before and to = a/to- In this case the corresponding potential V^'(x) 
has the same centrifugal behaviour at the origin. The analogous analysis of wave 
functions and of their transformations by symmetry operator can be performed. 

One can check that F(x) does not factorize the space dependence from a constant 
matrix. And also we stress that the potential V^'(x) cannot be made diagonal by 
global rotation and therefore it is not to be viewed as a pair of standard scalar 
reflect ionless potentials: this means that there is flux from one channel to the other. 

In order to ascertain the reducible or irreducible character of the solutions of fl6~4]) 
it is important to clarify the conditions for reducibility for the case = and 
VW(x) a transparent potential. As a consequence of (|3~9"D, fl4l"l) , (f42|) these conditions 
read: 

V® = W 2 -W = -AF'; (67) 
V® = W 2 + W' + A = 0; (68) 
W 2 + W = W 2 - W' + A; (69) 
W + W = 2F; (70) 
B = WW -W' = 2WF + 2F' (71) 



9o{x) = 
9i(x) = 

92(x) = 
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with solution 

A = -e; W = i- W = 2F-C 

providing as a consequence the condition [F, F'\ = 0, corresponding to the factor- 
ization of x— dependence of the matrix F(x). It is easy to check that our potentials 
V^(x) do not satisfy these conditions and therefore are not reducible. 

4. Conclusion 

We have demonstrated that the SUSY approach allows to relate as SUSY part- 
ners a dynamical matrix system where the symmetry is manifest to another matrix 
system where this symmetry is hidden in the sense that it is not otherwise easy to 
guess it (see e.g. Subsection 3.1.). Therefore one is connecting systems with more 
complex dynamics to systems which are simpler or even solvable. For Matrix Quan- 
tum Mechanics this approach provides examples of dynamical (matrix) systems and 
associated symmetry operators: it is useful since a straightforward general investi- 
gation of symmetries of a dynamical systems is not an easy task. In particular the 
connection between degeneracy of levels and the existence of symmetry operators 
is not mandatory and needs further clarification which presumably will depend on 
detailed dynamical properties of the system under investigation. 

This line of research is not academic but on the contrary should become a useful 
approach to investigate non-trivial quantum mechanical systems. We have restricted 
ourselves to 1-dim Matrix Quantum Mechanics: the algebraic methods we develop 
seem therefore to be specifically suited for this type of dynamical systems for which 
spatial symmetries (like 0(3)) already have been used to reduce the problem to a one 
dimensional (radial) problem (separation of variables). 

In absence of general theorems we have studied first and second order intertwin- 
ing relations between matrix 1-dim hamiltonians discussing reducible and irreducible 
transformations among them. The examples we have discussed can be interpreted 
as coupled channel problems or Pauli-type Hamiltonians and our techniques may be 
instrumental to their diagonalization. We have provided for the first time explicit ex- 
amples of irreducible transformations in the context of transparent matrix potentials 
in the framework of Matrix SUSY Quantum Mechanics. 
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